Draft version February 2, 2008 

Preprint typeset using L^T^X style emulateapj v. 14/09/00 



ANALYTIC LIGHTCURVES FOR PLANETARY TRANSIT SEARCHES 

Kaisey Mandel 1-2 and Eric Agol 1,3 



ABSTRACT 

We present exact analytic formulae for the eclipse of a star described by quadratic or nonlinear limb darkening. 
In the limit that the planet radius is less than a tenth of the stellar radius, we show that the exact lightcurve can 
be well approximated by assuming the region of the star blocked by the planet has constant surface brightness. 
We apply these results to the HST observations of HD 209458, showing that the ratio of the planetary to stellar 
radii is 0.1207 ±0.0003. These formulae give a fast and accurate means of computing lightcurves using limb- 
darkening coefficients from model atmospheres which should aid in the detection, simulation, and parameter 
fitting of planetary transits. 

Subject headings: eclipses — occultations — stars: binaries: eclipsing — stars: planetary systems 



1. INTRODUCTION 

The eclipse of the star HD 209458 by an orbiting planet was 
recently used to measure the size and mass of the planet, which 
had been found with velocity measurements (Charbonneau et 
al. 2000; Henry et al. 2000). With this landmark discovery, 
the planetary transit tool was added to the planet-finder's tool- 
box, already yielding several planetary candidates (Udalski et 
al. 2002a; 2002b; Dreizler et al. 2002). Several large surveys 
which aim to find planets using the transit signature are now 
being carried out or planned and will soon yield large num- 
bers of lightcurves requiring fast computation of eclipse mod- 
els to find the transit needles within the haystack of variabil- 
ity (Borucki et al. 2001; Howell et al. 2000; Mallen-Ornelas 
et al. 2002; Koch et al. 1998; Deeg et al. 2000; Street et al. 
2002). Lightcurve fits to transit events may be used to char- 
acterize the planet and star, yielding important constraints on 
planet formation (Cody & Sasselov 2002; Hubbard et al. 2001; 
Seager & Mallen-Ornelas 2002). The data require an accu- 
rate description of limb-darkening as demonstrated by Hubble 
Space Telescope observations of HD 209458 of such high qual- 
ity that a quadratic limb-darkening law was needed to fit the 
transit lightcurve rather than the usual linear limb-darkening 
law (Brown et al. 2001). The limb-darkening of main-sequence 
stars is represented by functions of /i = cos 9, where 9 is the an- 
gle between the normal to the stellar surface and the line of sight 
to the observer (Figure la). Claret (2000) has found that the 
most accurate limb-darkening functions are the quadratic law 
in n and the "nonlinear" law which is a Taylor series to fourth 
order in p, 1 ' 2 ; the latter conserves flux to better than 0.05%. 

In this paper, we compute analytic functions for transit 
lightcurves for the quadratic and nonlinear limb-darkening 
laws, and make available our codes for the community (§ 7). 
For treatment of subtler effects during planetary transits see 
Seager et al. (2000); Seager & Sasselov (2000); Hubbard et 
al. (2001); Hui & Seager (2002). In section 2 we review the 
lightcurve of a uniform spherical source. In section 3 we derive 
the lightcurve for eclipses of nonlinear limb-darkened stars. In 
section 4 we give a simpler form in the limit of a quadratic 
limb-darkening law. In section 5 we give an approximation for 
the lightcurve in the case p < 0.1 which is very fast to compute 
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and is fairly accurate. In section 6 we apply the results to some 
example cases, and in section 7 we conclude. 

2. UNIFORM SOURCE 

We model the transit as an eclipse of a spherical star by an 
opaque, dark sphere. In what follows, d is the center-to-center 
distance between the star and the planet, r p is the radius of the 
planet, r» is the stellar radius, z = d/r* is the normalized sepa- 
ration of the centers, and p = r p /r* is the size ratio (Figure lb). 
The flux relative to the unobscured flux is F. 

For a uniform source, the ratio of obscured to unobscured 
flux is F e (p,z) = 1 - X e (p,z) where 




A c (p,z) = 



-(l+z 2 -p 2 ) 2 



l+p<Z 

\l-p\ <z< l + p 



p 
1 



z< 1- 
z< p- 



(1) 

and Ki =cos- 1 [(l-p 2 + z 2 )/2z], n = cos -1 [(p 2 + z 2 - \)/2pz]. 
We next consider the effects of limb-darkening. 

3. NON-LINEAR LIMB DARKENING 

Limb darkening causes a star to be more centrally peaked in 
brightness compared to a uniform source. This leads to more 
significant dimming during eclipse and creates curvature in the 
trough. Thus, including limb-darkening is important for com- 
puting accurate eclipse lightcurves. Claret (2000) proposed a 
nonlinear limb-darkening law that fits well a wide range of stel- 
lar models and observational bands, I(r)= l-X)«=i c n(l _ A i "' /2 X 
where /i = cos 9 = VT-r 2 , < r < 1 is the normalized radial 
coordinate on the disk of the star and I{r) is the specific inten- 
sity as a function of r or p with 7(0) = 1. Figure 1(a) shows the 
geometry of lensing and the definition of /i. The lightcurve in 
the limb-darkened case is given by 

- 1 d[F<(p/r,z/r)^ 
drl(r) — - 



F(p,z) = 
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dr 
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where F e (p,z) is the lightcurve of a uniform source defined in 
the previous section. 
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FIG. 1. — (a) Geometry of limb-darkening. Star is seen edge on, with the 
observer off the top of the page. The star has radius r* and 9 is defined as 
the angle between the observer and the normal to the stellar surface, while 
p, = cos0. (b) Transit geometry from perspective of observer. 

In what follows cq = 1 — c\ — C2 — C3 — C4. For convenience, 
we define a = (z-p) 2 , b = (z + p) 2 , and fl = X^=o c "( n + 4) _1 - 
We partition the parameter space in z and p into the regions and 
cases listed in Table 1 . Next we describe each of these cases in 
turn. 

In Case I the star is unobscured, so Fj = 1 . In Case II the 
planet disk lies on the limb of the star, but does not cover the 
center of the stellar disk. We define 
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In the above equations, B(a,b) is the Beta function, 
Fi(a,bi,b2,c;x,y) is Appell's hypergeometric function of two 
variables, and 2Fi(a,b,c;x) is the Gauss hypergeometric func- 
tion. The relative flux is F= l-(27rfi) _1 Y? n=0 Nc n {n+4Y l . This 
case covers the ingress/egress where the lightcurve is steepest. 
For Cases III and IV, the planet's disk lies entirely inside the 
stellar disk, but does not cover the stellar center. We define 



M = (l-a) (n+4) / 4 
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and L = p 2 (l-p 2 /2-z 2 
F = 1 - (40)- 1 



Then the relative flux is given by 



This case 



c p 2 + 2J2ti Mc n(n + 4)" 1 + c A L 
requires the planet to be less than half of the size of the 
star. In Case V the edge of the planet touches the center of 
the stellar disk, and the planet lies entirely within the stel- 
lar disk. The relative flux is F v = 1/2 + (2f2) _1 X^=o c "(" + 
4) -1 2^i (l/2,-(« + 4)/4,l;4/? 2 ). For Case VI the planet's di- 
ameter equals the star's radius, and the edge of the planet's disk 
touches both the stellar center and the limb of the star. The 
relative flux is 



«=0 



(5) 



In Case VII the edge of the planet's disk touches the stellar cen- 
ter, but the planet is not entirely contained inside the area of the 
stellar disk. The relative flux is 
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In Case VIII the planet covers the center and limb of the 
stellar disk. The relative flux is F = -f2 -1 Y^t=o c nN(n + 4)" 1 . 
This and the previous case apply when the planet is larger than 
half the size of the star. For the special Case IX the planet's 
disk lies entirely inside the stellar disk, and the planet cov- 
ers the stellar center. The relative flux is F = (4f2) -1 [co(l - 
p 2 ) + C4(l/2-L)-2J2l =l Cn(n + 4T l M]. This is the bottom 
of the transit trough for nearly edge-on inclinations if p <C 1 . 
In Case X the planet is concentric with the disk of the star, 
at the precise bottom of the transit trough. In this case, F = 
fT 1 J2 4 „=Q c n( l -p 2 ) (n+4)/4 (n+4T l . This formula applies only for 
edge-on orbits when there is a central transit. Finally, in Case 
XI the planet completely eclipses the star, so that F = 0. In this 
case the "planet" is likely a star. In the event that c\ = C3 = 0, 
these lightcurves can be simplified as we describe in the next 
section. 

4. QUADRATIC LIMB DARKENING 

In this section we describe the limb-darkening with a func- 
tion which is quadratic in p,, I(r) = 1 -71(1 -^-72(1 -yii) 2 , 
where 71 +72 < 1 ■ The nonlinear law in the previous section re- 
duces to this case when c\ = C3 = 0, 02 = 71+272, and C4 = -72. 
In this limit, the hypergeometric functions reduce to elliptic in- 
tegrals which are much faster to compute, so in this section we 
provide these simpler formulae. 

For a quadratic limb-darkening law, the lightcurve is 
F =l-(4ny l [(l-C2)\ e + c 2 (\ d + 2 : Q(p-z))-c 4 r] d ], where 
\ e is defined in Equation 1, while \ d and rf are given in Table 
1. 

In Table 1 , the various functions are 
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m = ^(p 2 + 2z 2 ). 



(7) 



where k = y/ (1 - a) / (4-zp) and q = p 2 -z 2 . Here U(n 7 k) is the 
complete elliptic integral of the third kind with the sign con- 
vention of Gradshteyn and Ryzhik (1994). For linear limb- 
darkening, 72 = 0, Merrill (1950) presents an equivalent an- 
alytic expression in terms of an "eclipse function," a. The 
expressions here require fewer evaluations of the elliptic inte- 
grals, which decreases computation time, and include quadratic 
limb darkening. Our expression for eclipse with quadratic limb- 
darkening decreases computation time by more than an order of 
magnitude compared to evaluating the expressions in the previ- 
ous section or numerical integration of the unocculted flux. 




5. SMALL PLANETS 

For a small planet, p < 0.1, the interior of the lightcurve, z < 
l-p, can be approximated by assuming the surface-brightness 
of the star is constant under the disk of the planet, so that F = 



( ,, , and /*(z) = (4-zp) 1 f*_ p I(r)2rdr. If one knows the 
limb-darkening coefficients of the star in question (from, say, 
spectral information), and if the semi-major axis is much larger 
than the size of the star so that the orbit can be approximated by 
a straight line, then the shape of the eclipse for p < 0. 1 is simply 
determined by the smallest impact parameter, zo = a p cosi/r*, 
where a p is the semi-major axis and i is the inclination angle. 
For example, at the midpoint of the eclipse, z = Zo, while at the 
1/4 and 3/4 phases of the eclipse, z = Zi/ 4 = (1 + 3z^) 1/2 /2. Tak- 
ing the ratio of the depth of the eclipse at these points yields 
R = (1 -F(zo))/(l -F(zi /4 )) = /(zo)//(zi/ 4 ). This then deter- 
mines an equation for zo, and the resulting zo can then be used 
to determine p = y/4il(l -F(zo))/I(zo). Whenl-/?<z< l+p 
and p < 0.1 an approximation to the lightcurve is 
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where /*(z) = (1 -a) 1 j z _ p I(r)2rdr. which is accurate to better 

than 2% of 1-F(0) for p = 0.1 and ^ n=l c n < 1 (see Figure 2). 
This is a very fast means of computing transit lightcurves with 
reasonable accuracy, and may be used for any limb-darkening 
function, generalizing the approach of Deeg et al. (2001). 



6. DISCUSSION 

Figure 2 shows five lightcurves, the first of which has c„ = 
0, {n = 1,4}, while the other four have c n = 1 , c m = 0, {m n} 
for p = 0. 1 ; these may be thought of as a basis set for any nonlin- 
ear limb-darkening. Note that the higher order functions have 
flux which is concentrated more strongly toward the center of 
the star, and thus have a more gradual ingress and a deeper min- 
imum as more flux is blocked at the center than the edge. All 
of the curves cross near z ~ 0.7, which means that accurate 
observations are required near minimum and egress/ingress to 
constrain the coefficients of the various basis functions. If the 
inclination is large enough that z > 0.7 for the entire transit, 
then it may be difficult to constrain the c„'s. 



FIG. 2. — Transit lightcurves for p = 0.1 and c\ = c 2 = c 3 = c 4 = (solid 
line), and all coefficients equal zero but c\ = 1 (dotted line), C2 = 1 (dashed 
line), C3 = 1 (dash-dot line), or c 4 = 1 (dash, triple-dot line). The lighter lines 
(nearly indistinguishable) show the approximation of section 5. 

As Claret (2000) claims that the nonlinear limb-darkening 
law is the most accurate, we have compared the transit with 
p = Q.\ for the nonlinear and quadratic laws. Of the entire 
grid of models computed by Claret (2000) which covers 2000 
K <T eff < 50000 K, < logg < 5, -5 < [M/H]< 1, and filters 
u, b, v, y, U, B, V, R, I, J, H, K, the largest difference between 
the quadratic and nonlinear models for p = 0.1 is 3% of the 
maximum of 1 —F. Thus, the quadratic law should be sufficient 
for main-sequence stars when an accuracy of less than 3% is 
required; indeed, the average difference for the entire grid of 
models is about 1% of the maximum value of 1 — F for each 
lightcurve. In absolute terms, this is about lO^ip/O.l) 2 of the 
total flux, an accuracy which can be achieved from space. 

So far we have only presented the lightcurve as a function of 
z and p. To determine z as a function of time requires the plan- 
etary orbital parameters, which for zero eccentricity is given 

1 /2 

as z = a p r~ l [(smut) 2 + (cos i cos cot) 2 ] , where u is the orbital 
frequency, while t is the time measured from the center of the 
transit. Contribution of flux from the planetary companion or 
other companions may be added to the lightcurve, reducing the 
transit depth. 

To illustrate the utility of our formulae, we have fit the non- 
linear limb-darkened lightcurve to the HST-STIS data of HD 
209458 (Brown et al. 2001). The best-fit parameters in this 
case are p = 0.12070 ±0.00027, i = 86.591° ±0.055, a p /r* = 
8.779 ±0.032, c\ = 0.701, c 2 = 0.149, c 3 = 0.277, c 4 = -0.297, 
with a reduced \ 2 = 1-046. The c„'s are poorly constrained 
given the small differences in the basis functions relative to 
the observed errors and the large impact parameter for this sys- 
tem. The errors on the other parameters are marginalized over 
the c„'s. Limiting the limb-darkening to quadratic, we find 
7i = 0.296 ± 0.025 and 72 = 0.34 ± 0.04, consistent with the 
values derived by Brown et al. (2001) and with stellar atmo- 
sphere predictions. The value for p in the quadratic case is 
consistent with the nonlinear case, indicating that the fit is inde- 
pendent of the assumed limb-darkening law. For a stellar mass 
of 1.1 ±O.1M (Mazeh et al. 2000) and period of T = 3.5248 
days (Brown et al. 2001), we find r» = 1.145 ± 0.035fl Q and 
r p = 1.376 ±0.043/? 7Kp . 

We apply the small-planet approximation described in the 
previous section to HD 209458, assuming quadratic limb- 
darkening. Cody & Sasselov (2002) determined the effec- 
tive temperature and surface gravity for this star, which imply 
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71 = 0.292 and 72 = 0.35 for the ]-band flux (close to the ef- 
fective wavelength of the HST STIS data) from the models of 
Claret (2000). From the lightcurve, one finds F(zo) = 0.9835 
and F(z]/ 4 ) = 0.9847. Solving the equations from Section 5 
gives zq = 0.546 and p = 0.12, very similar to the parameters 
derived by a fit to the entire lightcurve. This technique may be 
used for finding initial parameters for lightcurve fitting. 

7. CONCLUSIONS 

We have derived analytic expressions for an eclipse including 
quadratic limb-darkening and nonlinear limb-darkening. The 
nonlinear law (Section 3) provides an accurate fit to realistic 
stellar limb-darkening, while the quadratic fit (Section 4) pro- 
vides a fast means of obtaining a relatively accurate lightcurve. 
For an extremely fast and fairly accurate approximation for any 
limb-darkening law, the equations in Section 5 may be used to 
derive lightcurves. If the limb-darkening law is known from 
the spectral type of the star, then one can use the formulae in 
Section 5 to analytically estimate both the minimum impact pa- 
rameter (in units of stellar radius) and the ratio of the planetary 
radius to the stellar radius. We have written a code which takes 
the properties of a host star, finds the limb-darkening coeffi- 
cients in the tables of Claret (2000), and computes lightcurves 
for the parameters of a given planetary transit. This code will 



be useful for simulating planetary transit searches (Gaudi 2000; 
Defay et al. 2001; Jenkins et al. 2002; Remund et al. 2002; 
Jenkins 2002; Pepper & Gould 2002), searching for planetary 
transit signals in lightcurves collected by a given search, and for 
fitting and measuring the errors of the parameters of detected 
planetary transit events. Planetary searches suffer from two 
important backgrounds: grazing eclipsing binaries and triple 
systems in which two stars eclipse while the flux from the 
third reduces the depth of the eclipse. Using the appropriate 
limb-darkening coefficients for each star's spectral type will 
help to distinguish these contaminants from true planetary tran- 
sits, which can be accomplished using the formulae presented 
here. These routines are made available by download from 
http://www.pha.jhu. edu/^agol/. 
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